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ABSTRACT 

We present a deep, narrow J-band search for proto-galactic Lya emission at 
redshifts z ~ 9. We combine an exceptionally deep image of the Hubble Deep Field 
South, obtained using a narrow band filter centred on the wavelength 1.187/im using 
the VLT/ISAAC facility, with existing deep, broad band images covering optical to 
near infrared wavelengths. We search for z ~ 9 Lya-emitting galaxies displaying a 
significant narrow band excess relative to the Js~band that are undetected at opti- 
cal wavelengths. We detect no sources consistent with this criterion to the 90% point 
source flux limit of the NB image, Fnb = 3.28 x 10~^^ ergs s^^ cni^^. The survey selec- 
tion function indicates that we have sampled a volume of approximately 340 Mpc'^ 
to a Lya emission luminosity of lQ*^h~^ ergss"^. We conclude by considering the po- 
tential implications for the physics of the high-redshift universe. 
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1.1 The high— redshift universe 

Studies of galaxies and QSOs at the highest known red- 
shifts have historically provided some of the strongest con- 
straints upon theories of galaxy formation and the physical 
state of the Inter-Galactic Medium (IGM) . Observations of 
the high-redshift universe have been extended progressively 
during the past decade via a series of dedicated studies: 
Steidel et al. (1996) demonstrated that optical colour se- 
lection techniques, tuned to identify the spectral signature 
of young star forming galaxies attenuated by intervening 
neutral hydrogen, are highly effective at identifying large 
numbers (hundreds per square degree) of galaxies at red- 
shifts 2.5 < z < 4.5. Spectroscopic observations of large 
samples (approximately 1000 to date) of so-called 'Lyman- 
Break' galaxies has permitted computation of the galaxy 
luminosity function (Steidel et al. 1999), composite stellar 
populations (Pettini et al. 2001; Shapley et al. 2001) and 
the amplitude of large scale structure variations (Giavalisco 
et al. 1998). Employing similar colour selection techniques, 
samples of distant galaxies have been extended to redshifts 
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z — 6, with the advent of the Advanced Camera for Sur- 
veys (ACS) facility deployed on the Hubble Space Telescope 
(HST; Dickinson et al. 2004; Stanway et al. 2004a). Though 
such distant surveys constrain the volume averaged star for- 
mation rate (SFR) for a flux limited sample at early cosmic 
times (Giavalisco et al. 2004) , detailed statistical conclusions 
await the compilation of larger samples. Although the de- 
tection of a Lya-emitting galaxy at a redshift z = 10 was 
recently reported by Pello et al. (2004) following a colour- 
selected survey for redshift z > 7 galaxies located behind 
well-studied gravitational lens clusters, the exact nature of 
the source remains contentious (e.g. Weatherley et al. 2004; 
Bremer et al. 2004). 

Observational programmes employing narrow band 
(NB) imaging techniques combine the advantage of wide 
area imaging surveys with enhanced sensitivity to discrete 
spectral features. Though NB selection techniques restrict 
such programmes to line-emitting galaxies at specific red- 
shifts, spectroscopic confirmation of sources pre-selected to 
to display strong spectral features is relatively straightfor- 
ward compared to continuum selected samples. A number 
of NB imaging programmes have successfully extended the 
study of high-redshift objects from redshifts z = 3 — 4 
(Cowie and Hu 1998), through redshifts 2 = 4-6 (Mal- 
hotra and Rhoads 2002; Rhoads et al. 2003; Hu et al. 2004) 
to redshifts z = 6.5 (Hu et al. 2002a; Hu et al. 2002b; Ko- 
daira et al. 2003). Such studies have extended current con- 
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straints upon the star formation history of the universe to 
redshifts z ^ 6, albeit incorporating a bias toward emission 
line sources (i.e. unobscured star formation). 

Studies of high-redshift QSOs (of greater apparent 
brightness than galaxies at the same luminosity distance 
yet presenting surface densities lower by several orders of 
magnitude to a given flux limit) have been revolutionised 
by the Sloan Digital Sky Survey (SDSS; Fan et al. 2003). 
Detailed study of the current sample of ^ 10 QSOs identi- 
fied at redshifts z > 5.7 have provided compelling evidence 
for the onset of neutrality of the hydrogen component of the 
IGM (Becker et al. 2001), the existence of sustained metal 
production in QSOs (Pcntcricci ct al. 2002) and the pres- 
ence of massive black holes in the universe approximately 
0.9 Gyr after the epoch of recombination (Fan et al. 2001). 

1.2 H I absorption: the Lya forest and the IGM 

Though evidence exists that star forming structures are 
present in the universe at redshifts z > 6, the possibility 
to observe such structures is dependent critically upon the 
opacity of intervening material. The observed spectral en- 
ergy distributions (SEDs) of high redshift sources display 
Lyman series absorption arising from two distinct phenom- 
ena: discrete structures and continuous, diffuse absorption. 
Discrete H I structures along the line-of-sight to background 
sources each give rise to an individual Lya absorption fea- 
ture. As the absorption redshift increases, the discrete 'for- 
est' of Lya absorption merges to create a blanketed region 
of near contiuous absorption at wavelengths blueward of 
the Lya transition in the background source. It is this ab- 
sorption that gives rise to the Lyman-Break discontinuity 
present in the SEDs of all high-redshift sources. Songaila 
and Cowie (2002) estimate the amplitude of the Lyman 
Break discontirmity to be 

Am = 3.8 + 20.3 logio (^) ' (1) 

from the average transmission of the Lya forest region in 
the spectra of a sample of 15 QSOs located over the red- 
shift interval 4.42 < z < 5.75. However, extrapolation of 
this relation to redshifts 2; > 6 is complicated as the op- 
tical depth of the intervening IGM potentially approaches 
and exceeds unity. The high redshift IGM can be approxi- 
mated to a continuous spatial distribution of hydrogren gas. 
A uniform IGM becomes opaque (r ^ 1) to wavelengths 
A < 1216A at a neutral hydrogen fraction of 10~^ (Gunn and 
Peterson 1965). As the fraction of neutral gas increases, the 
absorption saturates and develops a characteristic 'damp- 
ing wing' - extending the absorption profile to wavelengths 
rcdward of Lya and potentially attenuating Lya emission in 
the embedded source. Spectroscopic observations of a subset 
of high-redshift QSOs identified within the SDSS (Becker et 
al. 2001) indicate the presence of absorption troughs where 
the linc-of-sight transparency is zero, i.e. apparently indi- 
cating r ^ 1. However, associating extended regions of zero 
transmitted flux with the Gunn-Peterson (GP) effect is com- 
plicated by the requirement to disentangle the effects of line 
blanketing in the Lya forest from the excess absorption aris- 
ing from the onset of neutrality in the IGM itself (Songaila 
and Cowie 2002). Although the above observations may have 
detected initial indications of the onset of neutrality in the 



IGM, it is important to note that, due to density fluctua- 
tions in the IGM and clustering of sources of ionising ra- 
diation, the apparent onset of neutrality in the IGM may 
be a strong function of the line-of-sight. Further observa- 
tions arc required to determine the extent of variations in 
the properties of the IGM due to these sources of cosmic 
variance (White et al. 2003). Furthermore, the observation 
of a Lya emitting galaxy at redshift z = 6.56 (Hu et al. 
2002a) is not inconsistent with observation of a possible GP 
effect in lower redshift QSOs. Haiman (2002), Santos (2004) 
and Barton et al. (2004) describe the ionising effect of a 
star forming galaxy embedded in a neutral IGM and note 
that, depending upon the exact assumptions made regarding 
the mass and star formation properties of the source and the 
physical conditions present in the IGM, the galaxy will form 
a local H 11 region of sufficient size to permit transmission of 
a partially attenuated Lya line and associated continuum. 



In contrast to observations of absorption troughs in 
redshift a ~ 6 QSO spectra, observations of the Cosmic 
microwave Background (CMB) have the potential to con- 
strain the properties of the universe during the reionisation 
epoch optical depth regime r 3> 1. Observation of a positive 
correlation in the Wilkinson Microwave Anisotropy Probe 
(WMAP) CMB polarisation-temperature map at large an- 
gular scales can be interpreted as a scattering event with 
an optical depth t = 0.17 - consistent with the reionisa- 
tion of the universe at redshifts 2: < 30 (Kogut et al. 2003). 
Although the exact variation of the mean neutral hydrogen 
fraction as a function of redshift is strongly model depen- 
dent, the initial WMAP results support an early commence- 
ment of the reionisation epoch. 



The present paper describes a dedicated search for 
high-redshift star forming galaxies, employing an extension 
of broad and narrow-band selection techniques applied at 
optical wavelengths to the near infrared (NIR) wavelength 
regime. In particular we focus upon the application of a 
narrow J-band filter centred at A = 1.187n.m to detect the 
signature of Lya emitting galaxies located about a redshift 
z = 8.8 (termed z ^ 9 in the following text). The present 
work builds upon earlier searches for galaxies at extreme red- 
shift employing a similar narrow J-band approach (Parkes 
et al. 1994), in addition to narrow Jf-band studies poten- 
tially sensitive to galaxies at even greater redshift (Bunker et 
al. 1995). The following sections are organised as follows: in 
Section 2 we describe the observational strategy, the choice 
of target field and we describe the narrow J-band data ac- 
quisition and reduction. We further describe the archival 
data employed and the creation of a multi-colour catalogue 
for the field. In Section 3 we discuss source detection, com- 
pleteness and contamination within the multi-colour cata- 
logue. In Section 4 we describe candidate emission line ob- 
jects detected in the field and we present our conclusions is 
Section 5. Unless otherwise indicated, values of f^M.o = 0.3, 
i^A.o = 0.7 and Ho = 70 kms~^ Mpc^^ arc adopted for the 
present epoch cosmological parameters describing the evo- 
lution of a model Friedmann-Robertson-Walker universe. 
Where used, h is defined as h = Hq/IQQ kms~^ Mpc~^. 
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2 OBSERVATIONS 

2.1 Observing strategy and field selection 

Following the discussion in Section 1.2 it is clear that can- 
didate redshift z ^ 9 galaxies will display a well defined 
photometric signature in a search employing optical to NIR 
broad band photometry with a narrow band centred on Lya 
emission. The combination of Lya forest and IGM absorp- 
tion will result in effectively zero flux transmission at rest 
wavelengths A < 1216A. However, uncertainty regarding the 
physical state of the IGM and the potential for z ~ 9 sources 
to create local H ii regions ensures that a range of Lya emis- 
sion properties may result from a given star forming galaxy 
population. Throughout this paper, we consider that a red- 
shift 2: ~ 9 Lya-emitting galaxy will display a significant 
narrow band excess relative to the J^-band, in addition to 
displaying a continuum break consistent with almost com- 
plete attenuation of photons at rest frame A < 1216A. In 
order to generate an effective survey for such sources cer- 
tain additional factors must be considered: 

(i) NIR continuum imaging data must achieve a limiting 
depth of AB ^ 25.5 . The brightest emission line galax- 
ies confirmed at redshifts z = 5.7 display AB magnitudes^ 
z' ^ 24.5 - 25 (Hu et al. 2004)^. The additional distance 
modulus between a redshift z = 5.7 and z = 8.8 results in 
an additional dimming term of 1 magnitude. 

(ii) Optical imaging data must reach a limiting depth 
typically 1.5 magnitudes fainter than NIR data. Early-type 
galaxies located at redshifts z ^ 2 can generate a spectral 
discontinuity between optical and NIR continuum bands of 
amplitude D ~ 1.5 mag. (Stanway et al. 2004b). Failure 
to identify the continuum break directly could lead to the 
mis-identification of redshifted [OII]3727 emission in such 
sources as candidate z ~ 9 Lya emission. 

Following these considerations the Hubble Deep Field 
South (HDFS; Williams et al. 2000) Wide field Planetary 
Camera 2 (WFPC2) apex pointing (a = 22''32'"55!464, S = 
-60°33'05'.'01, J2000) was selected as the target field in 
order to exploit the high quality of optical to NIR image 
data available for the field. In particular, the combination 
of HDFS WFPC2 and Very Large Telescope (VLT) Infrared 
Spectrometer And Array Camera (ISAAC; Moorwood 1997) 
observations of the field provide images to typical depths 
AB=28 and AB=26 in optical and NIR bandpasses respec- 
tively (Labbe et al. 2003). 

2.2 Narrow— band near infrared observations 

Narrow J-band observations of the HDFS WFPC2 pointing 
were obtained during ESQ Period 69 (May 19th to Septem- 
ber 17th 2002) employing the ISAAC facility mounted on the 
Nasmyth-B focus of the 8.2 meter VLT Antu telescope. The 

^ Source magnitudes are computed using the AB system (Oke 
1974), i.e. AB = -48.6 - 2.51ogF^, where is the spec- 
tral energy density within a particular passband in units of 
ergss^^ cm^"^ Hz^^. 

^ Note that we employ this comparison as, at redshifts z = 5.7 
and z = 8.8, redshifted H I absorption lies at the blue edge of the 
z' and Js bandpasses respectively. 
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Figure 1. The spectral response of the ISAAC NB119 filter su- 
perposed upon the J-band emission spectrum of the night sky. 
The night sky emission spectrum is taken from Rousselot et al. 
(2000). 

ISAAC short wavelength camera is equipped with a Rock- 
well Hawaii 1024 x 1024 HgCdTe array. The pixel scale is 
0'.'1484 pix"^ and the field size is 2.'5 x 2'5. The read noise 
of the detector is 16.7 electrons and the gain is 4.5 elec- 
trons ADU~^. Observations were performed employing the 
NB119 filter (Figure 1). The filter has an effective width of 
89. 5A and is centred at a wavelength 1.187/ttm (correspond- 
ing to the location of Lya emission at a redshift z = 8.76). 
The NB119 filter is particularly effective for the detection of 
faint emission as it is located in a region of low sky noise, 
largely isolated from the 'forest' of bright OH emission that 
dominates the sky background in the J-band. 

The science programme was allocated 32 hours of ser- 
vice mode observations, to be obtained in seeing conditions 
of 0'.'6 or better. Observations were scheduled as 32 one 
hour observing blocks (OBs), each consisting of ten, 300 sec- 
ond exposures dithered spatially according to a Poissonian 
random offset located within a box size of 20" in right as- 
cension and declination. The choice of pointing centre and 
dither sequence were made to produce a final NB image 
of the HDFS that was well matched to the available deep, 
broad NIR band observations obtained as part of the Faint 
Infra-Red Extragalactic Survey (FIRES) project (Labbe et 
al. 2003). The final data set consists of 42 hours of narrow 
band observations. This represents all usable data, includ- 
ing observations obtained outwith the specified observing 
constrains (sky transparency, moon phase and atmospheric 
seeing). The data consists of 420 individual science images 
corresponding to a total exposure time of 126,000 seconds. 

2.2.1 Data reduction 

Narrow band imaging data were a) corrected for vary- 
ing pixel response using twilight sky exposures, b) sky- 
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subtracted having masked array regions containing objects 
detected above a specified ADU level, c) corrected for both 
high- and low-frequency spatial artefacts, d) shifted to a 
common pixel scale and coadded using a suitable pixel re- 
jection algorithm and bad pixel mask. The applied data re- 
duction techniques are broadly similar to those described 
in Labbe et al. (2003) for deep JsHKs observations of the 
HDFS and are summarised below. 

Individual narrow band exposures were dark subtracted 
using the master dark frame obtained during each observing 
night. The flat field response and bad pixel masks were con- 
structed using the flat routine in ECLIPSE^ . Twilight flat 
field exposures were obtained employing a fixed exposure 
time and display uniformly decreasing/increasing count lev- 
els throughout a given sequence. Exposures displaying count 
levels less than 5000 ADUs were excluded as computation 
of pixel responses based upon these frames display a charac- 
teristic discontirmity at the location of the array read junc- 
tion indicating the presence of a residual dark signature. The 
residual dark signature represents an additive effect that can 
be minimised by constructing the fiat field frame employing 
twilight exposures displaying ADU levels greater than 5000. 
The pixel response is defined as the gradient of a linear re- 
gression fit to the set of values returned by a given pixel over 
the twilight sequence (the gain) normalised by the median 
gain of the array. Pixels displaying gain values outwith the 
interval [0.5,2] were defined as deviant pixels and added to 
the bad pixel mask. Flat field frames constructed from indi- 
vidual twilight sequences obtained over ESO Period 69 were 
compared. The pixel-to-pixel root mean square (rms) devi- 
ation between flat field frames is at the level of 0.2 — 0.4% 
per pixel. Large scale gradients in individual frames do not 
exceed 2%. A master flat field frame was then constructed 
from all suitable twilight exposures and employed to correct 
the pixel response in the science exposures. 

The NB119 filter samples an exceptionally low sky back- 
ground and relatively long individual exposure times are re- 
quired to ensure that resulting images are limited by Poisso- 
nian sky noise and not by detector read noise. The median 
sky value in a typical 300 second exposure obtained with 
the NB119 filter is 300 ADU pixel"^ where the exact value 
varies most significantly as a function of airmass (Figure 
2). Large variations in total sky emission or the presence of 
gradients are not noted within a typical one hour observing 
sequence (10 science frames). Accurate sky subtraction can 
therefore be achieved employing standard techniques em- 
ploying a measure of the sky in a given frame drawn from a 
statistical combination of science frames obtained adjacent 
in time. The sky contribution in each pixel within a given 
science frame is computed from the six adjacent frames in 
time. The sky contribution per pixel is estimated via scaling 
the six frames to a common median level and computing the 
average pixel value after rejection of the highest and lowest 
values. 

Pixel values in the sky image created for each frame 
are biased to higher ADU levels by the contribution of as- 
tronomical sources. Therefore, the location of bright objects 
within individual images must be identified with an object 



^ ECLIPSE is an image processing package written by N. Devil- 
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Figure 2. Median sky value determined for 420 science frames as 
a function of airmass computed at the mid-point of each frame. 
The error bars represent the root mean square deviation of the 
sky pixel distribution in each frame. 



mask and removed from second computation of the sky level. 
Relative integer pixel shifts generated by the dithered ob- 
serving sequence are computed using the relative location 
of bright stars in the sky subtracted images. The 420 indi- 
vidual images are registered and co-added (averaged) having 
rejected the 10 highest and lowest pixels to produce an inter- 
mediate image. All astronomical sources displaying a peak 
ADU level greater than a 5 cr threshold are identified within 
the CO added image and a set of masklets are created cor- 
responding to the location of all objects in the dithered im- 
ages. The above sky subtraction procedure is then repeated 
excluding pixels identified by the object masklet associated 
with a particular frame. 

Sky subtracted images contain low level residuals (of the 
order of a few ADUs) at both high and low spatial frequen- 
cies. High frequency residuals arise from a residual dark sig- 
nature (also termed pattern noise), the amplitude of which is 
related to gradients in the dark signal, i.e. it is largest at low 
row numbers (1,2,...) and those following the read junction 
(513, 514, ...). This signature was removed by computing the 
median value of all pixels contributing to a given row, having 
first masked both bad pixel and object locations. Analysis 
of corrected images indicated a residual signature associated 
with the columns of each image. Therefore, the residual dark 
signature in each column was corrected for via application of 
the same procedures. Low frequency residuals largely arise 
from imperfect sky subtraction and are removed by fitting a 
slowly varying function (2nd order cubic spline) along rows 
and columns incorporating bad pixel and object masking. 

At this stage, background corrected, sky subtracted im- 
ages contain cosmic ray events and additional bad pixels not 
identified within the flat field computation. Cosmic rays are 
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identified employing the PHIIRS^ routine crzap: each image 
is median smoothed over a given spatial scale and the result 
is subtracted from the original image. Application of a sigma 
threshold based upon the sky noise in each frame identifies 
the signature of cosmic ray events. Events thus identified 
are added to the bad pixel mask associated with that frame. 
Unidentified bad pixels remaining in each frame are identi- 
fied via analysis of the rms-normalised pixel distribution in 
each frame with bad pixel and object masking applied. Pix- 
els displaying absolute rrns normalised values greater than 
six are added to the bad pixel mask associated with the 
particular frame. This procedure typically identifies 40 — 50 
additional bad pixels per ten frame observing sequence. Im- 
ages analysed in this manner were then registered employing 
integer pixel shifts and co-added employing both a bad pixel 
mask and a pixel rejection criteria that rejected the 10 high- 
est and lowest values associated with a particular pixel. A 
pixel rejection criteria was still found to be necessary in or- 
der to exclude deviant pixels associated with objects (11% 
of the image area) from the co-added image. 

The quality of the coadded image was improved em- 
ploying a weighting scheme to optimize the image combina- 
tion given variations in seeing, sky transparency and back- 
ground noise within individual images. Individual shifted 
images, d, were combined to produce a final image, F, em- 
ploying a weighting scheme such that F — '^Wid, where 
Wi = Zi/{visl). The applied weighting scheme respectively 
corrects for relative variations in sky transparency via Zi, 
derived from the mean 5" diameter aperture flux of four 
bright reference stars within the HDFS field, variations in 
background noise via, Vi, the measured variance per pixel 
and variations in atmospheric seeing, Si, measured from the 
profiles of the same four bright reference stars. The distri- 
bution of image weights are normalised such that ^ = 1. 
The Pull-Width at Half-Maximum (FWHM) of the Point 
Spread Function (PSF) in the reduced image is 0'.'47. The 
dimensions of the final coadded image arc 2. '9 x 2.' 9. A sci- 
ence image was created from this frame based upon pixels 
displaying weighted exposure times greater than a fraction 
0.95 of the 126,000 on-sky exposure time. The dimensions 
of this science image are 2'.2 x 2'.2. 



2.2.2 Flux calibrating the NB119 image 

The photometric zero point for the NB119 image was com- 
puted employing 47 bright, isolated sources common to the 
NB119 and FIRES Ja-band images. Assuming that the 
SEDs of these calibration sources display no strong disconti- 
nuities, the Js-band fiux density provides an accurate esti- 
mate of the flux density with the NB119 filter, i.e. Js — NB = 
0. Source detection was performed on the NB119 image 
and apertures were transformed to the FIRES Js image 
to compute ,Js magnitude values. The transformation was 
computed using the IRAF^ routines geomap and geoxytran. 



Pat Hall's Infrared Imaging Reduction Software. 
® IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which axe operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 3. Computation of the zero— point of the NB119 image. 
The dotted line indicates a linear relationship of the form y = 
X + ZP. Individual error bars are displayed but possess a mean 
amplitude 0.006 magnitudes. 



Twenty bright, rmclcatod galaxies were employed as refer- 
ence sources and the root mean square deviation of the ref- 
erence sources about the computed transformation was 0.2 
pixels in the x- and j/-pixel directions. Source brightness 
measures wore computed within 5" diameter circular aper- 
tures and the distribution of Js versus NB119 instrumental 
magnitude was found to be linear with a gradient of unity. 
The relationship is displayed in Figure 3 and indicates that 
1 ADU received in the NB119 image corresponds to an 
AB magnitude of 22.64 ± 0.02. 



2.3 Archival optical and near infrared 
observations 

Optical to NIR broad band photometry of the HDFS is taken 
from the catalogue of Labbe et al. (2003). The authors em- 
ploy as input material for the multi-colour catalogue seven 
images, corresponding to the U3ooB45oVeo6l8i4JsHKs fil- 
ters, that have each been transformed and resampled to the 
astrometric system defined by the F814W HST image and 
have been convolved with a kernel that reproduces the seeing 
of the FIRES -ff-band image as determined via the median 
FWHM of bright stars (FWHM = 0'.'48.). Photometric mea- 
sures are computed employing 0'.'7 diameter apertures and 
are zero pointed on the AB magnitude system. As noted in 
Section 2.2.2, the NB119 image is not transformed to the 
astrometric system defined by the FIRES catalogue images 
as this would require a resampling of the data and the un- 
necessary introduction of correlated noise between adjacent 
pixels within our detection image. Instead, the set of 0'.'7 di- 
ameter detection apertures generated by the NB119 image 
(see Section 3 below) are transformed to the FIRES astro- 
metric system employing the transformation described in 
Section 2.2.2. Finally, no additional kernel is applied to cor- 
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rect for the slight diflterence in the measured seeing between 
the NB119 and FIRES //-band image as the measured dif- 
ference (AFWHM = Q/.'Ol) is smaller than the 0'.'04 variation 
in seeing over the //-band image reported by Labbe et al. 
(2003). 



3 SOURCE DETECTION AND 

COMPLETENESS 

Source detection and analysis was performed on the NB119 
image using the SExtractor software package (Bertin and 
Arnouts 1996). Circular apertures of diameter 0'.'7 were em- 
ployed to compute source fluxes. Source detection criteria 
were tuned to select objects displaying at least 12 pixels 
above a threshold of OAffsky Prior to source detection, the 
image was convolved with a 'mexican hat' filter of FWHM 
0'.'55. As a band-pass filter, the mexican hat emphasizes 
sources of extent comparable to its width, i.e. close to the 
PSF extent in this case. This was deemed an acceptable 
choice given the assumption that target a ~ 9 sources will 
appear unresolved in both the NB119 and Js images. 

The flux completeness limit of the NB119 image was 
estimated by introducing and recovering artificial sources 
within the field. The same analysis was performed upon 
the FIRES Js band image for comparison. Rcdshift z ^ 
9 sources arc assumed to be unresolved in both images. 
The simulation proceeded by sub-dividing each image into 
28x28 cells of 5" side length. Each of the 784 cells was fur- 
ther divided into 256 regularly spaced points and an artificial 
source was constructed at each point in the grid. All simu- 
lated sources in a particular image are introduced with the 
same source profile and flux. The image containing the arti- 
ficial sources is then added to the original image. Simulated 
sources possess low brightness values (AB > 24) and their 
contribution to the total photon noise per pixel is negligible. 
Source extraction is then performed on the image containing 
simulated sources. To avoid image blending among the simu- 
lated sources, the source simulation and image addition plus 
recovery operation is performed for each of the 256 locations 
within the grid of 5" cells in the simulated frame. For each 
5" cell the detection probability is the ratio of recovered to 
introduced sources. 

Variations in the detection probability within a given 
map are dominated by the local efi^ect of bright objects 
within the image on the likelihood to detect and extract sim- 
ulated sources. The simulation of 256 source positions within 
each 5" cell generates a grid of source centroids separated 
by O'.'S. The presence of a bright object within a particular 
5" cell will effectively mask all sources simulated in that cell 
and reduce the corresponding detection probability. The ef- 
fect is highly localised and results in the apparently noisy 
aspect in the detection probability maps. These detection 
probability variations are of considerable importance to any 
measure of the clustering of faint sources in the field but do 
not compromise the estimation of the mean completeness 
within a given source fiux interval. The mean value of the 
detection probability at a given fiux level, computed over 
the entire field, corresponds to the completeness as a func- 
tion of source fiux and is displayed for the NB119 and Js 
images in Figure 4. Adopting the 90% point source recovery 
threshold as the limiting magnitude in each band generates 




AB magnitude 

Figure 4. Mean detection probability of simulated point sources 
as a function of AB magnitude within the NB119 (square sym- 
bols) and Js images (triangular symbols). The horizontal dashed 
line indicates the 90% completeness threshold. 



magnitude limits of NB ^ 25.2 and Js ^ 26.2 respectively. 
The rms amplitude of sky counts measured in the NB image 
within the applied 0'.'7 aperture corresponds to a magnitude 
AB = 28.05. This figure is approximately 20% higher than 
the rms variation anticipated from the extrapolation of back- 
ground counts in raw data frames. The integrated signal- 
to-noisc ratio (SNR) of an NB detection at AB = 25.2 is 
therefore 13.8. The corresponding SNR value of a J^-band 
detection at AB = 26.2 is approximately 10. The detectabil- 
ity of such apparently significant sources is primarily af- 
fected by the additional photon noise contributed by bright 
sources within the image area. We note that the computed 
90% point source completeness limit of Js = 26.2 agrees 
well with the corresponding value of Js = 26.3 derived by 
Labbe ct al. (2003) given the slight difl'erences between the 
appUed source simulation method and SExtractor source 
detection parameters in each case. In common with Labbe 
ct al. (2003) we present 0'.'7 aperture photometry throughout 
this paper. The correction required to convert photometric 
measures computed in 0'.'7 apertures to 5" apertures (which 
we assume to be 'total' measures) was determined to be 0.7 
magnitudes via analysis of bright stars in the NB image. 

In addition to determining the completeness properties 
of the NB119 image, the false detection rate was also esti- 
mated. A set of 100 artificial images, containing no objects, 
and displaying random background noise properties consis- 
tent with the final NB119 image, was constructed. Source 
extraction was performed for each artificial image and the 
mean false detection rate was computed. With the combina- 
tion of SExtractor parameters described above, a false de- 
tection rate of 0.7 events per NB119 'noise' image resulted. 
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Figure 5. Distribution of NB excess Js — N3 versus NB mag- 
nitude for all objects in the NIR field area of the HDFS (solid 
points). The solid horizontal line indicates Js — ATS = 0. The ver- 
tical dotted line indicates the selection criterion NB ^ 25.2. The 
horizontal dotted line indicates the selection criterion — NB ^ 
0.3 and corresponds to a Lya EW of 21 K (Appendix A). The 
dashed curve indicates the predicted J a — NB uncertainty (3 a) 
as a function of NB magnitude returned by the completeness 
analysis (Section 3). 



4 RESULTS 

Figure 5 displays NIR narrow band excess, expressed as Js — 
NB, versus NB magnitude for all sources extracted from 
the total area covered by the optical and NIR data for the 
HDFS field (Figure 6). Sources displaying - NB > 0.3 
and NB < 25.2 are flagged as potential 2 ~ 9 emitting 
galaxies. This NB excess figure corresponds to a rest frame 
Lya equivalent width (EW) of approximately SA (Appendix 
A) and represents a 3 cr deviation from the flat spectrum 
hypothesis at the NB magnitude threshold {NB < 25.2) - 
itself the 90% point source completeness limit (Figure 4). 

Two sources are identified by the above selection cri- 
teria, HDFS-352 and HDFS-768, where the identification 
numbers refer to the HDFS photometric catalogue of Labbe 
et al. (2003). The photometric spectrum formed by the 
U3ooB4,5oVfjOf>Isi4JsHKs plus NB photometry of each object 
is displayed in Figure 7. In each case the addition of deep 
optical photometry is sufficient to exclude each NB excess 
detection as a potential « ~ 9 source: each NB excess object 
is detected in all optical bands and the photometric rcdshift 
of each source (Rudnick et al. 2001) - Zphot ~ 1-54 ± 0.06 



(HDFS-352) and Zphot = 0.76! 



(HDFS-768) 



sistent with a NB excess arising from redshifted H/J4861 
and Ha 6563 emission respectively. Although no candidate 
2^9 sources arc detected within the HDFS survey area, the 
confirmation of two sources whose photometric redshifts are 
consistent with the observed NB excess arising from narrow 
line emission lends support to the assertion that the adopted 



Figure 7. Low resolution spectral energy distributions formed by 
the multi-band photometry of the two NB excess objects identi- 
fied within the HDFS. The clear detection of each source within 
optical passbands precludes each source as a candidate z ~ 9 
source. 



selection criteria identify faint, narrow emission line galax- 
ies. 

The Lya emission line selection function generated by 
the NB excess search technique is described in Appendix 
A. Although the applied NB < 25.2 magnitude threshold 
corresponds to a total flux integrated across the NB fil- 
ter of Fnb = 3.28 X 10^^* ergs s~^ cm~^, computation of 
the flux transmitted from a Lya-emitting galaxy located 
at a redshift z and displaying a rest frame velocity width 
ffv, requires a more detailed investigation of the applied 
selection criteria, the spectral response of the NB119 fil- 
ter and the assumed Lya source properties. Figure A2 dis- 
plays the co-moving volume sampled as a function of Lya 
emission luminosity for three values of the rest frame veloc- 
ity width of the Lya-emitting source. For the case where 
the rest frame velocity width of putative z ~ 9 sources is 
cr„ = 50kms~^, the NB excess survey area samples a co- 
moving volume of 340 Mpc^ to a Lya emission luminos- 



ity of LLya ^ 10" 



' h ^ ergs s ^ . 



5 CONCLUSIONS 

No candidate Lya emitting sources at redshifts 2 ~ 9 have 
been identified within the HDFS survey area by the selection 
criteria described in Section 4. Expressing the observational 
selection criteria in terms of the volume of the universe sam- 
pled to a given Lya luminosity permits an investigation of 
this null result in terms of the instrinic space density and 
luminosity of putative Lya emitting sources at z 9. It 
must be noted however, that we currently only consider the 
case where the Lya emission from « ~ 9 sources is unob- 
scured, i.e. the escape fraction of Lya photons is unity. This 
scenario corresponds to the case where the Lya absorption 
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Figure 6. Greyscale NB119 image of the HDFS. The image orientation is North up and East left. The field geometry of the NIR broad 
band data set is identical to the NB image. The field geometry of the HST WFPC2 field is indicated. Note that the Planetary Caraera 
(shaded region) data does not contribute to the final catalogue. The total field area contributing to the NB excess catalogue is 4 arcmin^ . 
The location of the two NB excess sources identified from the catalogue are indicated. 



opacity due to the intervening IGM, and to the Lya emitting 
source itself, is zero. Decreasing the relative fraction of Lya 
photons that escape from z ~ 9 sources will increase the 
corresponding upper limit on the unobscured volume aver- 
aged luminosity density derived from the NB excess survey. 
However, the degeneracy that exists between the absorption 
properties of the IGM and the unobscured versus observed 
volume averaged Lya luminosity density at z ~ 9 compli- 
cates a more complete investigation of the null result gen- 
erated by the current NB excess survey. The range of star 
formation conditions required to maintain an ionised IGM 
at a redshift a ~ 9 remains a subject of active discussion (c.f. 
Haiman 2002; Santos 2004) and is not considered further in 
this paper. 

A more empirical approach is to consider known Lya- 
emitting sources and to transpose their properties to redshift 
z ^ 9. Hu et al. 2004 (hereafter H04) present a search for 
Lya emitting galaxies at 2 ~ 5.7 employing a narrow band 
filter centred on the wavelength 8150A. Though the H04 
study does not present the highest redshift Lya sources cur- 
rently known (c.f. Kodaira et al. 2003), their study is based 
upon near complete spectroscopic follow up observations - 
19 out of 23 observed sources (from a total sample of 26 can- 
didates) were confirmed as Lya at z = 5.7 - and the overall 
design of their study is very close to the NIR selected survey 



presented in this paper. H04 report a surface density of Lya 
emitting galaxies of 0.03 arcmin"^ to a Lya flux limit of 
2 X 10~^^ ergs s~^ cm~^ from a total areal coverage of 700 
arcmin"^ . 

Within our adopted cosmological model, an unresolved 
source at a redshift z = 8.8 is 2.7 times fainter than an 
identical source viewed z = 5.7. A further factor of 0.5 must 
be applied to account for the different fraction of total flux 
measured by the 3" diameter apertures applied by H04 and 
the 0'.'7 diameter apertures employed in the current study. 
The H04 flux limit therefore corresponds to an approximate 
limit of 3.5 x 10~^* ergs s~^ cm~^ when transposed to a 
redshift z = 8.8 and employing 0'.'7 diameter apertures, i.e., 
assuming that no additional evolution occurs, the z = 5.7 
sources observed by H04 are sufficiently bright to be ob- 
served within our NIR selected survey. Applying the sur- 
face density of confirmed z = 5.7 Lya emitting galaxies to 
the areal coverage of the current NIR survey, indicates a 
probability to detect & z ^ 9 galaxy of 0.12 - assuming no 
evolution between the two redshifts. This 'volume shortfall' 
indicates that the current NIR selected survey will have to 
be extended by up to eight times the currently sampled area 
in order to realistically probe the no evolution scenario. 

The present narrow J-band selected survey has suc- 
cessfully demonstrated that the sensitivity required to de- 
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tect putative z ^ 9 Lya emitting galaxies emission has been 
achieved for the case where no evolution occurs between the 
redshifts z = 5.7 and z = 8.8. However, it is clear that, in 
order to place stronger hmits on the space density of such 
sources, the arcal coverage of the current study must be ex- 
tended. These observations are currently in progress and will 
be presented in a future paper. 
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APPENDIX A: THE EMISSION LINE 
SELECTION FUNCTION FOR THE NB119 
FILTER AND THE VOLUME SAMPLED AS A 
FUNCTION OF Lya EMISSION LUMINOSITY 

The survey sensitivity to Lya emission at a redshift z ~ 9 
may be quantified by considering the volume sampled as a 
function of Lya emission luminosity. In order to constrain 
the volume averaged Lya emission at redshifts 2: ~ 9, a de- 
tailed understanding of the emission line selection function 
generated by the NB119 filter is required, i.e. the probabil- 
ity to identify an emission line of given total flux, velocity 
width and central wavelength. 

The flux ratio implied by the Js — NB ^ 0.3 selec- 
tion criterion corresponds to an oberved frame Lya equiv- 
alent width limit of or 2.8A in the rest frame of a 
z = 8.8 source. For comparison, the observed frame EW 
limit employed by H04 is approximately IIOA. The lower 
limit achieved in the current study results from a lower 
narrow- minus reference-band selection limit and the use 
of a narrower filter. The applied NB ^ 25.2 magnitude 
threshold corresponds to a total flux integrated across the 
NB filter of Fnb = 3.28 x 10"^* ergss"^ cm'^. The above 
EW threshold therefore indicates that the Lya emission flux 
contributing to this NB flux is approximately one-third of 
the total received flux, i.e. fhya = 1.1 x 10^^** ergss^^ cm^'^. 
The flux sensitivity calculation is applicable to continuum 
sources of spectral slope similar to the mean slope of bright 
calibration sources within the HDF South (stars and bright 
galaxies), i.e. Js — NB — (Section 2.2.2). However, com- 
putation of the emission line selection function requires the 
colour term describing the fractional difference in flux re- 
ceived from an emission line of particular velocity width 
and central wavelength (assuming a Gaussian line proflle) 
compared to the moan calibration source. The NB119 filter 
transmits approximately 85% of the total flux received from 
a narrow emission line located at the central wavelength of 
the fllter. The fraction of the incident flux received from a 
Gaussian emission line of varying velocity width and central 
wavelength, T(Ac,(t(,), may be computed by convolving the 
emission line profile with the spectral response of the NB119 
filter, i.e. 



T(Ac 



J^F{\-K,al)S{\)d\ 
/„°°F(A-Ac,ai)dA ' 



(Al) 



where F{X — Ac,o-^) is a Gaussian lino profile of central 
wavelength, Ac, and observed-frame velocity width, a'y. The 
NB119 fllter response is given by S{\). The quantities of 
central wavelength and velocity width are readily expressed 
as redshift, z, and rest frame velocity width, assuming 
that the emission is Lya. Contours describing r(z,CT„) are 
displayed in Figure Al. We do not consider the case where 
the Lya emission proflle is non-Gaussian. 

As noted above, a source satisfying J — NB ^ 0.3 and 



The effective widtfi of the NB119 filter is 85. SA. We do not 
consider tlie attenuation of the putative x ~ 9 emission spec- 
trum by the intervening IGM in this calculation. However, as the 
NB119 spectral response lies very close to the blue wavelength 
cut-oflF of the Js filter, partial absorption of the Js flux compared 
to the NB119 flux is not considered to be a major systematic 
uncertainty. 
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Figure Al. Fraction of flux transmitted, T(z, ay), for a Gaussian 
Lya emission profile of redshift, z, and rest frame velocity width, 
(Tv, convolved with the spectral response of the NB119 filter. The 
horizontal dotted line indicates the redshift at which Lya emission 
lies at the central wavelength of the NB119 filter. 



NB ^ 25.2 corresponds to a transmitted Lya emission flux 
fhya ^ 1.1 X lO"'^* ergss"'^ cm~^. The incident Lya flux, 
fiya ) is equal to this transmitted flux divided by the trans- 
mitted fraction, i.e. 



fhya — f'Lya/T{z,av). 



(A2) 



The flux limit for a line of fixed rest frame velocity width 
corresponds to a varying Lya emission luminosity limit as a 
function of redshift within the assumed cosmological model 
according to Li,ya{z) = Att dl{z)flya{z,(7v), where dLiz) is 
the luminosity distance as a function of redshift. The co- 
moving volume sampled as a function of Lya emission lumi- 
nosity, V(LLya) is then expressed as 
/• 00 

V(LLyc.) = 6dn / H(LLya)(d\/(z)/dz)d«, (A3) 
Jq 

where b corrects for the photometric completeness (6 = 0.9), 
df2 is the survey solid angle, and 



H(LLya) = 1 where LLya ^ LLya(x,0-„) 

= otherwise. 



(A4) 



Figure A2 displays the volume sampled as a function of Lya 
emission luminosity by the NB119 image presented in this 
paper. Three rest frame velocity width values arc considered 
to demonstrate the dependence of the luminosity versus vol- 
ume sampled values upon the assumed velocity width. The 
analysis indicates that the NB survey selection criteria sam- 
ple a volume of approximately 340 h~'^ Mpc"^ to a Lya emis- 
sion luminosity of 10'*^ ergs in the case where the 
rest frame velocity width of the source is ffy = 50 kms"'^. 

This paper has been typeset from a T^^X/ file prepared 

by the author. 
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Figure A2. Co— moving volume sampled as a fmiction of Lya 
emission luminosity. Three values of the rest frame Lya velocity 
width are displayed; <t„ = 20 (solid line), 50 (dashed line) and 
100 (dotted Une) kms"!. 



